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A nearby core collapse supernova will produce a burst of neutrinos in several detectors worldwide. 
With reasonably high probability, the Earth will shadow the neutrino flux in one or more detectors. 
In such a case, for allowed oscillation parameter scenarios, the observed neutrino energy spectrum 
will bear the signature of oscillations in Earth matter. Because the frequency of the oscillations in 
energy depends on the pathlength traveled by the neutrinos in the Earth, an observed spectrum 
contains also information about the direction to the supernova. We explore here the possibility of 
constraining the supernova location using matter oscillation patterns observed in a detector. Good 
energy resolution (typical of scintillator detectors), well known oscillation parameters, and optimisti- 
cally large (but conceivable) statistics are required. Pointing by this method can be significantly 
improved using multiple detectors located around the globe. Although it is not competitive with 
neutrino-electron elastic scattering-based pointing with water Cherenkov detectors, the technique 
could still be useful. 

PACS numbers: 14.60. Pq, 95.55.Vj, 97.60.Bw 



I. INTRODUCTION 

The core collapse of a massive star leads to emis- 
sion of a short, intense burst of neutrinos of all fla- 
vors. The time scale is tens of seconds and the 
neutrino energies are in the range of a few tens of 
MeV. Several detectors worldwide, both current and 
planned for the near future, are sensitive to a core 
collapse burst within the Milky Way or slightly be- 
yond d]. 

The first electromagnetic radiation is not expected 
to emerge from the star for hours, or perhaps even 
a few days. Therefore any directional information 
that can be extracted from the neutrino signal will 
be advantageous to astronomers who can use such 
information to initiate a search for the visible su- 
pernova. We note that not every core collapse may 
produce a bright supernova: some supernovae may 
be obscured, and some core collapses may produce 
no supernova at all, in which case directional infor- 
mation will aid the search for a remnant. 

The possibility of using the neutrinos themselves 
to point back to the supernova has been explored 
in the literature [21 [3]. Triangulation based on rel- 
ative timing of neutrino burst signals was also con- 
sidered in [2.; however available statistics, as well 
as considerable practical difficulties in prompt shar- 
ing of information, makes time triangulation more 
diflicult. Leaving aside the possibility of a TeV 
neutrino signal [3 (which would likely be delayed), 
the most promising way of using the neutrinos to 
point to a supernova is via neutrino-electron elastic 
scattering: neutrinos interacting with atomic elec- 
trons scatter their targets within a cone of about 



25° with respect to the supernova direction. The 
quality of pointing goes as ~ N~^/'^, where N is 
the number of elastic scattering events. In wa- 
ter and scintillator detectors, neutrino-electron elas- 
tic scattering represents only a few percent of the 
total signal, which is dominated by inverse beta 
decay Ve + P n + e'^ , for which anisotropy is 
weak [31. Furthermore the directional information 
in the elastic scattering signal is available only for 
water Cherenkov detectors, for which direction in- 
formation is preserved via the Cherenkov cone of 
the scattered electrons. Taking into account the 
near-isotropic background of non-elastic scattering 
events [3], a Super- K-like detector |5] (22.5 kton fidu- 
cial volume) will have 68% (90%) C.L. pointing of 
about 6° (8°) for a 10 kpc supernova; this could im- 
prove to < 1° for next-generation Mton-scale water 
detectors. Long string water detectors [6] do not 
reconstruct supernova neutrinos event-by-cvcnt and 
so cannot use this channel for pointing. Scintillation 
light is nearly isotropic and so scintillation detec- 
tors have very poor directional capability, although 
there is potentially information in the relative posi- 
tions of the inverse beta decay positron and neutron 
vertices [3S], and some novel scintillator directional 
techniques are under development [7]. 

We consider here a new possibility: detectors with 
sufficiently good energy resolution will be able to 
obtain directional information by observing the ef- 
fects of neutrino oscillation on the energy spectrum 
of the observed neutrinos, assuming that oscilla- 
tion parameters are such that matter oscillations are 
present. Although not competitive with elastic scat- 
tering, some directional information can be obtained 
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even in a single detector (unlike for time triangula- 
tion). Combinations of detectors at different loca- 
tions around the globe may yield fairly high quality 
information. 



II. DETERMINING THE DIRECTION 
WITH EARTH MATTER EFFECTS 

Supernova neutrinos traversing the Earth's mat- 
ter before reaching a detector will experience matter- 
induced oscillations, depending on the values of the 
MNS matrix parameters [51 [TUl ITTl [T?l IT5]. 
Whether or not there will be an Earth matter effect 
depends on currently-unknown mixing parameters, 
013 and the mass hierarchy: matter oscillation will 
occur for both v^. and v^, for values of sin^ Oiy, <, 10^^, 
for normal but not inverted hierarchy; if ^13 is rel- 
atively large, sin'^ ^ 10"'^, then matter oscillation 
occurs for but not i^g for either hierarchy [14, 15]. 
The frequency of the oscillation in L/E, where E 
is the neutrino energy and L is the neutrino path- 
length in Earth matter, depends on now fairly well- 
known mixing parameters. Therefore, the oscillation 
pattern in neutrino energy E measured at a single 
detector contains information about the pathlength 
L traveled through the Earth matter. If the path- 
length L is known, one knows that supernova is lo- 
cated somewhere on a ring on the sky corresponding 
to this pathlength. If another pathlength is mea- 
sured at a different location on the globe, the loca- 
tion can be further constrained to the intersection of 
the allowed regions. 

A Fourier transform of the inverse-energy distri- 
bution [5] of the observed neutrinos will yield a peak 
if oscillations are present. References [5] and |^ ex- 
plore the conditions under which peaks are observ- 
able with a view to obtaining information about the 
oscillation parameters. The authors assume that the 
direction of the supernova, and hence the pathlength 
through the Earth, is known. Here we turn the ar- 
gument around: we assume that oscillation param- 
eters are such that the matter effects do occur and 
can be identified, and that enough is known about 
MNS parameters to extract information about L and 
hence about supernova direction from the data. A 
similar idea to determine possible georeactor loca- 
tion from the oscillated spectrum was explored in 
reference [T6|. We note that by the time a nearby 
supernova happens, the hierarchy and whether ^13 
is large or small may in fact be known from long- 
baseline and reactor experiments. With reasonably 
high probability [17], the Earth will shadow the su- 
pernova in at least one detector. We note that lack 
of observation of a matter peak in the inverse-energy 
transform (assuming there should be one) gives some 



direction information as well: if no peak is present, 
one can infer that the supernova is overhead at a 
given location. If the hierarchy and value of 6*13 are 
already known with sufficient precision at the time of 
the supernova, we should know in advance whether 
or not a peak in the k distribution should appear; 
otherwise, its appearance for at least one detector 
location may answer the question. 



III. EVALUATION OF THE CONCEPT IN 
IDEALIZED SCENARIOS 

To evaluate the general feasibility of this con- 
cept we make several simplifying assumptions. We 
consider only inverse beta decay in large water 
Cherenkov and liquid scintillator detectors (we ig- 
nore the presence of other interactions, which should 
be a small correction; some of them can be tagged) 
[55] . We will first consider a detector with per- 
fect energy resolution, and then consider resolutions 
more typical of real water Cherenkov and scintillator 
detectors. 

We borrow some of the assumptions and notation 
of reference [8]. We assume a neutrino interaction 
cross-section proportional to E^, perfect detection 
efficiency above threshold and no background. We 
assume a "pinched" neutrino spectrum of the form: 

average neutrino energy. We choose parameters 
a — 3, average energies for the flavors = 15 MeV 
and Ei,^ = 18 MeV, and |^ = 0.8. These param- 
eters correspond to the "Garching" model [18]. We 
ignore for this study "spectral sphts" {e.g. [T^ ) 
or other features which will introduce additional 
Fourier components. We assume that there are no 
non-standard neutrino interactions or other exotic 
effects that modify the spectra. 

The oscillation probabilities have been computed 
by numerical solution of the matter oscillation equa- 
tions [20^ using these vacuum parameters and the full 
PREM Earth density model [5T]. Between neigh- 
boring radial points in the model the matter den- 
sity is taken to be constant such that the three- 
neutrino transition amplitude may be computed fol- 
lowing the methods outlined in [25]. The final am- 
plitude is the product of all amplitudes across the 
matter slices along the neutrino's trajectory. The 
initial flux of neutrinos is taken to arrive at the 
Earth as pure mass states such that the detection 
probability is taken according to the probability of 
a neutrino being Pg flavor when it reaches the de- 
tector. The oscillation parameters were chosen to 
be sin^ 26*12 = 0.87, sin^ 26*13 = 0, sin^ 26*23 = 1-0, 
Am?2 = 7.6x 10-5 eV^ and Am^g = 2.4x10-3 eV^ 
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A. Perfect Energy Resolution 



The spectrum of inverse beta decay events, in- 
tegrated over time, is shown in Fig. [T] Fig. [T] 
shows on the bottom the "inverse-energy" spectrum, 
where the inverse-energy parameter y is defined as 
y = i2.5^cv ^ Fig. [2] shows the Earth matter mod- 
ulation of the spectrum, for L — 6,000 km. Shown 
on the bottom is the modulation in inverse-energy, 
for which the peaks are evenly spaced. 
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Figure 1: Top: assumed neutrino event spectrum with- 
out oscillations. Bottom: inverse-energy distribution. 

The Fourier transform of the detected inverse- 
energy spectrum is g{k) = f{y)e^'''^dy. The 
power spectrum GaF{k) — \g{k)\'^ assuming per- 
fect energy resolution is shown in Fig. [3] for no 
matter oscillation on the top and for matter os- 
cillation on the bottom, assuming pathlength L = 
6, 000 km. The power spectra are generated from the 
normalized inverse-energy distributions for which 
aF{y)dy = 1. Thus the power spectra are nor- 
malized so that G(jf(0) = 1. Fig. ijshows the power 
spectra for several values of L, illustrating how the 
peak moves to higher k values as the pathlength in- 
creases. For pathlengths such that the neutrinos tra- 
verse the Earth core {L > 10,700 km), additional 
peaks are present in the spectrum 0. There is no 





Figure 2: Top: assumed neutrino event spectrum with 
matter oscillations for L = 6, 000 km. Bottom: inverse- 
energy distribution with matter oscillations. 



observable peak for L less than about 2500 km, for 
which the neutrinos are no longer traversing much 
high-density matter. 

Fig. |5] shows now the effect of finite statistics, for 
a simulated supernova with 10,000 events and one 
with 60,000 events. The finite statistics result in 
a background for the main peak(s) in the power 
spectrum. For most of the following, we consider 
a rather optimistically large (but not unthinkable) 
60,000 event signal, which would correspond to a su- 
pernova at a distance of about 5 kpc observed with 
a 50 kton detector. 



1. Method for Determining Directional Information 

If one measures fcpoak, the position of the largest 
peak in the power spectrum, for a supernova signal, 
one can in principle determine the pathlength trav- 
eled by the neutrino in the Earth. We use a simple 
Neyman construction method f23] to estimate the 
quality of directional information. 

We first find the position of the largest peak in k 
as a function of pathlength L, assuming perfect en- 
ergy resolution but finite statistics. To find the peak 
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Figure 3: Inverse-energy power spectrum without (top) 
and with (bottom) matter oscillations. 




Figure 4: Examples of inverse-energy power spectra for 
several pathlengths. 



in the power spectrum, we first set a lower threshold 
of fc = 40 and an upper threshold of fc — 210. Below 
that threshold, the peak merges with the low k peak 
(corresponding to the unoscillated spectrum) and 
can no longer be identified. Peaks beyond k = 210 
would correspond to distances greater than the di- 
ameter of the Earth. For each k within that range we 
then evaluate the integral from k — Ak/2 to k+Ak/2 
which corresponds to the area under the peak. We 



Figure 5: Examples of inverse-energy power spectra for 
perfect energy resolution but finite statistics. Both lines 
show L = 6000 km. 



take the k for which this value is highest as the 
peak in the spectrum. We chose Afc = 4. Even 
though Fig. |4] suggests that peaks can be wider than 
that, we found more fluctuation in the peak's po- 
sition for higher Ak when taking finite energy res- 
olution into account, especially for small distances 
{L < 4000 km). Fig. [6] shows that the value of kp^ak 
is clearly correlated with pathlength L; for distances 
less than about 2000 km, for which the neutrinos do 
not undergo matter oscillations, it represents mainly 
random noise. The multiple peak structure for neu- 
trinos passing through the core is clearly visible for 
L > 10,700 km. We note that the height of the 
largest peak also contains information about L, as 
do the secondary peak positions, if such exist. 
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Figure 6: Distribution of the position of the maximum 
peak in as a function of matter-traversed pathlength 
L, assuming perfect energy resolution. There are 500,000 
simulated supernovae per L, each with 60,000 events. 

Given a particular measurement of fcpcaki one can 
then determine a range of distances L allowed, mak- 
ing use of the Neyman construction shown in Fig. [7] 
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To ensure contiguous regions in k we drop regions 
that contribute less than 3% to the final Neyman 
construction and increase existing regions instead so 
that the total covered area is 68% or 90%. The range 
in L values can then be mapped to an allowed re- 
gion on the sky. We have checked explicitly that the 
statistical coverage is as expected. 
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Figure 7: Neyman construction for fepeak and L: for a 
given measured fcpeak one reads off a range of allowed L 
values. The green area shows the 68% confidence region 
and the red area the 90% one. 



Fig. |8] shows an example Hammer projection sky 
map in equatorial coordinates showing 90% CL. al- 
lowed regions for an assumed true supernova direc- 
tion (indicated by a star) of R.A.= 20^ and decl.= 
-60° (occurring at 0:00 GMST), for assumed perfect 
energy resolution and statistics of 60,000 events. 

Fig.|9]shows the distribution of fractional sky cov- 
erage for perfect energy resolution. The distribution 
is bimodal, because the L < 2500 km possibility 
(corresponding to large fractional sky coverage) is 
often not excluded at 90% in the Neyman construc- 



tion. Fig. 10 shows the average sky coverage vs. dec- 
lination of the supernova, averaged over 24h of right 
ascension, for a detector located in Finland (63.66° 
N, 26.04° E). 

If we incorporate also information about the 
height of the largest peak h into a Neyman construc- 
tion, for long pathlengths we can remove the possi- 
bility of a short-pathlength overhead supernova, and 
improve the pointing quality significantly. Fig. [TT] 
shows the correlation between peak heights and L. 
Figs, [s] (bottom) and 10 show the effect of incorpo- 
rating this information. Subsequent plots will as- 
sume use of both Fourier peak position and height 
information. 




Figure 8: Example Hammer projection skymaps in equa- 
torial coordinates, showing 90% C. L. allowed regions on 
the sky for a supernova at the position indicated by a 
star. A 60,000 neutrino event signal measured in Fin- 
land with perfect energy resolution is assumed. The top 
plot shows the allowed region without taking into ac- 
count peak height; the bottom plot takes into account 
peak height information. 
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Figure 9: Histogram of fractional sky coverages for the 
90% CL. region, assuming perfect energy resolution in 
a single detector, using fcpeak information only, for the 
example configuration of Fig. [S] 



2. Combining Detectors 

Clearly, having several detectors around the globe 
observing the neutrino burst will improve the mea- 
surement. If each of the detectors could select a sin- 
gle L, an observation with two detectors will produce 
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better. 
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Figure 10: Sky coverage averaged over right ascension as 
a function of declination, for a single detector. The black 
line just takes into account the peak position, whereas 
the red line includes also the height of the peak. In total, 
83,500 supernovae, evenly distributed over declination, 
have been simulated for both lines. 
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Figure 11: Distribution of the heights of the maximum 
peak in fe as a function of matter-traversed pathlength 
L, assuming perfect energy resolution. There are 500,000 
simulated supernovae per L, each with 60,000 events. 



two allowed regions where the rings on the sky over- 
lap, and a third observation will narrow it down to 
one spot. However because more than one L region 
may be allowed for a given detector, the combination 
can include multiple regions. 

For the multiple detector case, we make the 
Neyman construction for 100,000 randomly chosen 
(k\\ hi, k2', /i2, ■ • .)-tuples only (with 100 bins in fcpeak 
and height for each detector) in order to compute 
it in a reasonable amount of time. In this case the 
smoothing procedure described above is not applied. 



Fig. 12 shows example scenarios involving two and 



three detectors and Fig. 13 summarizes average qual- 
ity as a function of declination. Clearly in this ideal- 
ized situation, combined information is quite good, 
and the more detectors spread around the globe, the 
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Figure 12: Combined skymaps for detectors with perfect 
resolution. Top: Two detectors with 60,000 events each. 
Bottom: Three detectors with 60,000 events each. 




Figure 13: Sky coverage averaged over right ascension 
as a function of declination, for one, two and three de- 
tectors. For one detector 83,500 supernovae have been 
simulated, for two detectors the number is 2,630 and for 
three detectors it is 760. 



B. More Realistic Detectors 

Next we will assume a slightly more realistic sit- 
uation. Imperfect energy resolution will tend to 
smear out the oscillation pattern and degrade the 



detectability of the peak in k. We estimate the ef- 
fect of energy resolution by selecting events from the 
spectrum and smearing their energies according to 
a Gaussian of the prescribed width. The energy res- 
olution functions used, the same as in reference 
are shown in Fig. |14[ one is characteristic of scin- 
tillator and one of water Cherenkov detectors. For 
water Cherenkov we assume a threshold of 5 MeV 
and for scintillator we assume a threshold of 1 MeV. 




10 20 30 40 50 60 7( 

Energy [MeV] 



Figure 14: Energy resolution functions used for water 
and scintillator. AE is the standard deviation cr of a 
Gaussian. 



1. Water Cherenkov Detectors 
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Figure 16: Distribution of the position of the maximum 
peak in as a function of matter-traversed pathlength 
L, assuming scintillator energy resolution. There are 
5,000,000 simulated supernovae per L, each with 60,000 
events. 



Clearly the water Cherenkov resolution smears the 
power spectrum information enough to preclude its 
use for this purpose; furthermore, far superior direc- 
tion information will come from elastic scattering in 
a water Cherenkov detector. Therefore we will focus 
subsequent attention on scintillator detectors, which 
have significantly better energy resolution and weak 
intrinsic direct pointing capabilities. 



Fig. 15 shows the distribution of fcpoak and L 
for simulated supernovae in a detector with water- 
Cherenkov-like energy resolution. Fig. [16] shows the 
same for scintillator. 
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Figure 15: Distribution of the position of the maximum 
peak in as a function of matter-traversed pathlength L, 
assuming water Cherenkov energy resolution. There are 
500,000 simulated supernovae per L, each with 60,000 
events. 



2. Scintillator Detectors 



Existing and near-future scintillator detectors 
with supernova neutrino detection capabilities are 
KamLAND [24], LVD [23 [26], Borexino [27] and 
SNO-I- [35]; these are however probably too small 
to acquire the large statistics required for this tech- 
nique. Future scintillator detectors of the tens of 
kton scale for which this technique could be feasi- 
ble are LENA [52], to be sited in Finland, and the 
ocean-based HanoHano [ 50] . 

Fig. [it] shows an example skymap for a scintillator 
detector located in Finland. Fig. [18] shows average 
sky coverage vs. declination for three examples of 
event statistics. 

Next we consider the case when multiple detectors 
are operating: Fig. [19] shows the results of combin- 
ing the information from two and three scintillator 
detectors located in Finland, ofl^ the coast of Hawaii 
(19.72° N, 156.32° W) and South Dakota (44.45° N, 
103.75° W). Fig. [20] shows average sky coverage vs. 
declination for these configurations. 
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Figure 17: Example scintillator skymaps, for a single de- 
tector located in Finland, assuming a 60,000 event signal 
(top). 
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Figure 18: Average scintillator sky coverage vs declina- 
tion for a single detector located in Finland, for 20,000 
event, 40,000 and 60,000 event signals. In total 200,000 
supernovae have been simulated for the 60,000 events 
case and 12,000 and 15,000 for the 20,000 and 40,000 
event cases, respectively. 



3. Incorporating Relative Timing Information 

We consider briefly now the possibility of incor- 
porating relative timing information between detec- 
tors to break degeneracies in the allowed region(s). 
A detailed study of the triangulation capabilities for 
specific neutrino signal and detector models is be- 
yond the scope of this work. We instead do some 
back-of-the-envelope estimates based on those in ref- 
erence [2J. For a signal registered in two detectors, 
the supernova direction can be constrained to a ring 
on the sky at angle 9 with respect to the line be- 
tween the detectors, with cos 9 = At/d and width 



(5(cos( 



5(At) 



where d is the distance between the 



detectors and d{At) is the time shift uncertainty be- 
tween the pulses. We assume S{At) ^ 30 ms/^/Ni, 
where A^i is ^1% of the total signal. A sharp feature 
in the signal timing could reduce 6{At). 

Fig. [21] shows an example for two detectors 
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Figure 19: Example scintillator skymaps, for a two (top) 
and three (bottom) detectors setup, each with a 60,000 
events signal per detector. 
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Figure 20: Average scintillator sky coverage vs declina- 
tion for one, two and three detectors with 60,000 events 
each. For one detector 200,000 supernovae have been 
simulated, for two detectors the number is 3,500 and for 
three detectors it is 1,256. 



(located in Finland and Hawaii), with the time- 
triangulated allowed region superimposed: the inter- 
section clearly narrows down the allowed directions. 

We can imagine also that another, non- 
scintillator, neutrino detector (or even a gravita- 
tional wave detector, e.g. j31j) could provide rela- 
tive timing information as well. For example, Ice- 
Cube at the South Pole could yield few ms iYoi- 
ing [T" ' 



Fig. 22 shows an example of the intersection 
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Figure 21: Example two scintillator detector skymap, 
with estimate of allowed region based on relative timing 
information superimposed (dark band). 



of the estimated IceCube plus single scintillator de- 
tector time-triangulation allowed region (assuming 
S{At) ~ 1 ms) with the single scintillator oscillation 
pattern region. 
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Figure 23: Effect of varying Ami2 on the peak position 
as a function of L; in each case other oscillation param- 
eters are held at their nominal values. 
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Figure 22: Example single scintillator detector skymap, 
with estimated allowed region determined from relative 
timing with the IceCube signal (dark band). 



IV. DISCUSSION 

We have assumed in these idealized scenarios per- 
fect knowledge of oscillation parameters. In prac- 
tice, imperfect knowledge of the oscillation param- 
eters will create some uncertainties. In particular, 
the power spectrum peak position is sensitive to the 
value of Amf2; the peak height is sensitive to both 
Aml^ and 6*12; 6*13 also has an effect on both 

"-peak 

and h. (The oscillation pattern is quite insensitive to 
the 23 mixing parameters.) Figs. 23 and 24 show the 
effect on fcpcak and h values of varying the oscillation 
parameters within currently allowed ranges [33J. 

From these plots one can infer that ^1% knowl- 
edge of the mixing parameters is desirable. However, 
one can be quite optimistic that such precision will 
have been attained by the time a core collapse su- 



Figure 24: Effect of varying 612 and 6^13 on the peak 
height as a function of L; in each case the other oscilla- 
tion parameters are held at their nominal values. 



pernova happens when a large scintillator detector 
is running. 

Another uncertainty that will affect the quality 
of pointing is that of the density of matter in the 
Earth. We found only small differences in kp^ak and 
h from varying the mantle density by ±3%, or from 
varying the overall density by ±5%, but observed 
some changes in peak pattern for the case of neutri- 
nos passing through the core when varying the core 
density by ±10%. 

Many other effects may degrade the quality of di- 
rection information that can be obtained using this 
technique. There may be real spectral features {e.g. 
"splits" ) which introduce additional Fourier compo- 
nents that could mask the peak, and detector im- 
perfections may do the same. We acknowledge also 
that there may be practical difficulties with the rapid 
exchange of information between experimenters re- 
quired for prompt extraction of directional infor- 
mation from multiple detectors. Nevertheless this 
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technique represents an interesting possibility- even 
half the sky is better than no directional infor- 
mation. The oscillation pattern gives information 
about direction with even a single detector, and en- 
hances any multiple-detector time-triangulation in- 
formation. Even if information from only a single 
detector is available, or if there are significant am- 
biguities, one can imagine also looking at the inter- 
section of the allowed region with the Galactic plane 
regions for which supernovae are most likely to occur 
([17]) (perhaps using the known probability distribu- 
tion as a Bayesian prior) to improve the chances of 
finding the supernova: see Fig. [25] 
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Figure 25: Average scintillator sky coverage vs declina- 
tion for a single detector in Finland, with the expected 
probability for supernova occurrence superimposed from 
reference [17| . 

We note that these estimates of pointing quality 
have been done using a fairly simple technique based 
on only two parameters characterizing the power 
spectra. One can imagine employing more sophis- 
ticated algorithms, e.g. making use of secondary 
peaks or matching to a template, and possibly in- 
corporating knowledge of specific detector proper- 
ties or neutrino flux spectral features. So although 
real conditions may degrade quality, with this sim- 
plified study we have not fully exploited all poten- 
tially available information. 



As a final note: the technique could in principle 
work to determine directional information for neu- 
trino signals from other astrophysical sources, such 
as black hole-neutron star mergers [34', assuming 
suflicient statistics. 



V. SUMMARY 

We have explored a technique by which experi- 
ments with good energy resolution can determine 
information about the direction of a supernova via 
measurement of the matter oscillation pattern. This 
method will only work for favorable (but currently 
allowed) oscillation parameters; it requires large 
statistics, good energy resolution, and well-known 
oscillation parameters, and it works best for rela- 
tively long neutrino pathlengths through the Earth. 
The method is especially promising for scintillator 
detectors. The criteria will be fulfilled in optimistic 
but not inconceivable scenarios. Combining infor- 
mation from multiple detectors, and possibly incor- 
porating relative timing information, may provide 
significant improvement. The method is inferior to 
that using elastic scattering in imaging Cherenkov 
(or argon time projection chamber) detectors; elas- 
tic scattering remains the best bet for pointing to 
the supernova. However it is possible that a super- 
nova will occur when no such detector is running, 
in which case one should use whatever directional 
information can be extracted from the observed sig- 
nals. 
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